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Abstract 

Several uniform asymptotics expansions of the Weber parabolic cylinder functions are 
considered, one group in terms of elementary functions, another group in terms of Airy 
functions. Starting point for the discussion are asymptotic expansions given earlier by 
F.W.J. Olver. Some of his results are modified to improve the asymptotic properties and 
to enlarge the intervals for using the expansions in numerical algorithms. Olver's results 
are obtained from the differential equation of the parabolic cylinder functions; we mention 
how modified expansions can be obtained from integral representations. Numerical tests 
are given for three expansions in terms of elementary functions. In this paper only real 
values of the parameters will be considered. 
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1 Introduction 

The solutions of the differential equation 

fy ( i 

dz 2 



\z 2 + a)y = U (1.1) 



are associated with the parabolic cylinder in harmonic analysis; see [20]. The solutions are 
called parabolic cylinder functions and are entire functions of z. Many properties in connection 
with physical applications are given in 

As in [|]], Chapter 19, we denote two standard solutions of (|l.l|) by U(a, z), V(a, z). These 
solutions are given by the representations 
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and the confluent hyper geometric function is defined by 

y> (a)n Z n 
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with (a) n = T(o + n)/r(a), n = 0, 1, 2, . . .. 
Another notation found in the literature is 

£>„(*) =Z7(-i/-|, z). 
There is a relation with the Hermite polynomials. We have 



hi z 
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2-«/ 2 e 3 ^„(z/V2), 
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F (n + J, z) = 2-™/ 2 e 3 " (-i) n H n (iz/^ 

Other special cases are error functions and Fresnel integrals. 
The Wronskian relation between U(a,z) and V(a,z) reads: 

U(a, z)V'(a, z) - U'{a, z)V(a, z) 



(1.2) 



(1.3) 



(1.4) 



(1.5) 



(1.6) 



which shows that U(a,z) and V(a,z) are independent solutions of ( |1.1| ) for all values of a. 
Other relations are 



U(a,z) 



V(a,z) 



TT 



cos 2 7rar(| + a] 

ra + a) 



[V(a, —z) — sin ira V(a, z)] , 



(1.7) 



7T 



[sin7ra U(a, z) + U(a, —z)] ■ 



The functions y\(a,z) and y 2 (a, z) are the simplest even and odd solutions of ( |1 . 1| ) and the 
Wronskian of this pair equals 1. From a numerical point of view, the pair {yi,y 2 } is not a 
satisfactory pair ||, because they have almost the same asymptotic behavior at infinity. The 
behavior of U(a, z) and V(a, z) is, for large positive z and z ^> \a\: 



U(a,z) 



1 7 2 

e z 



1 _2 



[l + 0(z- 2 )] 



V(a,z) =^2j^e l z 2 [l + O (z~ 2 )] . 



(1.8) 
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Clearly, numerical computations of U{a,z) that are based on the representations in (1.2) 
should be done with great care, because of the loss of accuracy if z becomes large. 

Equation Ql.ip has two turning points at ±2y/—a. For real parameters they become 
important if a is negative, and the asymptotic behavior of the solutions of (LI) as a —> — oo 
changes significantly if z crosses the turning points. At these points Airy functions are 
needed. By changing the parameters it is not difficult to verify that U(— %fi 2 , fity/2) and 

V(—^fi 2 ,fit\/2) satisfy the simple equation 



y ..4 (,2 



[t 2 - l) y = 0, (1.9) 



with turning points at t = ±1. For physical applications, negative a— values are most impor- 
tant (with special case the real Hermite polynomials, see ( |1.5|) ). For positive a we can use 
the notation U^/j?, jity/2) and V{\[i 2 , /xtv2), which satisfy the equation 

|f - {? + 1) y = o. (i.io) 

The purpose of this paper is to give several asymptotic expansions of U(a, z) and V(a, z) 
that can be used for computing these functions for the case that at least one of the real 
parameters is large. In pH] an extensive collection of asymptotic expansions for the parabolic 



cylinder functions as \a\ — > oo has been derived from the differential equation (1.1). The 
expansions are valid for complex values of the parameters and are given in terms of elemen- 
tary functions and Airy functions. In Section 2 we mention several expansions in terms of 
elementary functions derived by Olver and modify some his results in order to improve the 
asymptotic properties of the expansions, to enlarge the intervals for using the expansions in 
numerical algorithms, and to get new recursion relations for the coefficients of the expansions. 
In Section 3 we give similar results for expansions in terms of Airy functions. In Section 4 we 
give information on how to obtain the modified results by using integral representations of 
the parabolic cylinder functions. Finally we give numerical tests for three expansions in terms 
of elementary functions, with a few number of terms in the expansions. Only real parameters 
are considered in this paper. 

1.1 Recent literature on numerical algorithms 



Recent papers on numerical algorithms for the parabolic cylinder functions are given in [14] 
(Fortran; U(n,x) for natural n and positive x) and |l]| (Fortran; U(a,x),V(a,x), a integer 
and half-integer and x > 0). The methods are based on backward and forward recursion. 

H gives programs in C for U(a, x), V(a, x), and uses representations in terms of the con- 
fluent hyper geometric functions and asymptotic expressions, including those involving Airy 



functions. [23] gives Fortran programs for computing U(a, z),V(a, z) with real orders and 
real argument, and for half- integer order and complex argument, The methods are based on 
recursions, Maclaurin series and asymptotic expansions. They refer also to H for the evalu- 

ation of U(—ia, ze ) for real a and z (this function is a solution of the differential equation 
y" + (jz 2 ~ a )y = 0). [19] uses series expansions and numerical quadrature; Fortran and C 
programs are given, and Mathematica cells to make graphical and numerical objects. 

Maple ||] has algorithms for hypergeometric functions, which can be used in (IL^) and 



1.3Q . Mathematica [21] refers for the parabolic cylinder functions to their programs for the 



hypergeometric functions and the same advice is given in [12]. 
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2 Expansions in terms of elementary functions 

2.1 The case a < 0, z > 2</ = a, -a + 21 > 0. 

(Diver's expansions in terms of elementary functions are all based on the expansion 0-(4.3) [] 

V / (t 2 - 1)4 s=0 V 

as // — > 00, uniformly with respect to t G [1 + £,00); e is a small positive number and £ is 
given by 



e = §iv / ^T- I hi |t + V* 2 -lJ. (2.2) 

The expansion is valid for complex parameters in large domains of the fx— and t— planes; 
details on these domains are not given here. 

The coefficients A s (t) are given by the recursion relation 

1 1 dAJt) if* 3u 2 + 2 , , N , , ,. , . 

2 Vt 2 - 1 at * J Cs+1 (y? -1)2 

where the constants c s can be chosen in combination with the choice of g(fj,). Olver chose the 
constants such that 

= ^2 _ l)3s/2' ( 2 ' 4 ) 

where the u s (t) are polynomials in i of degree 3s, (s odd), 3s — 2 (s even, s > 2). The first 
few are 

, , tit 2 - 6) -9t 4 + 249t 2 + 145 
u (t) = l, Ul (t)= \ A , Mt) = — , 

and they satisfy the recurrence relation 

(t 2 - l)u' s (t) - 3stu s (t) = r s _i(i), (2.5) 

where 

8r s (i) = (3t 2 + 2)u s (t) - 12(s + l)tr,_i(i) + 4(t 2 - ly'^t). 
The quantity g(fJ-) in ( |2.1| ) is only available in the form of an asymptotic expansion: 

> (2-6) 

where 

1„2 1 1 ,,2 12 1 
^) = 2 ? " ? e^/ 3 , (2.7) 

50 = 1, 31 = ^, °3 = - 207360 ' 92s = (s = 1, 2, . . .), 

and in general 

o s = lim (2.8) 



1 We refer to Olver's equations by writing 0-(4.3), and so on. 



2.1.1 Modified expansions 

We modify the expansion in ( |2.lD by writing 



U(-^^tV2)=^^F,(t), F,(t)~f^, (2.9) 



where h(fx) and £ are as before, and 

t 



1 



(2.10) 



The coefficients 4> s (r) are polynomials in r, ^>o( r ) = 1) an d are given by the recursion 

</> 8 +i(t) = -4r 2 (r + l) 2 -^Mt) ~\f (20u 2 + 20u + 3) cf> s (u) du. (2.11) 

This recursion follows from (2.3) by substituting t = (r + |) j \Jt{t + 1), which is the inverse 
of the relation in (|2.1C| ). Explicitly, 

Mr) =1, 

01 (r) =-^(20r 2 + 30r + 9), 

02 ( T ) = (6160t 4 + 18480r 3 + 19404t 2 + 8028r + 945), (2.12) 

288 

r 3 

Mt) = -■fttttttC 27227200 ' 7 " 6 + l 22 522400r 5 + 220540320r 4 + 
5184L) 

200166120t 3 + 94064328r 2 + 20545650r + 1403325), 

where r is given in ( |2.10| ). Observe that limt_> 00 r(i) = and that all shown coefficients <j) s (j) 
vanish at infinity for s > 0. These properties of 4> s (t) follow by taking different constants c s 
than Olver did in (12.31). In fact we have the relation 



OO OO i f \ OO /.\ 

V V ^ S ^ T ' ~ V Us ^ > 

2-^, ,,2s 2-^, ,,2s ~ 2^, 



where the first series appears in ( |2.6| ), Explicitly, 

3 s 



(t) = (f-iy s Y J 9s- j <i> j (T)- ( 2 -i3) 

j=0 



The relation ( 2.13| ) can easily be verified for the early coefficients, but it holds because of the 
unicity of Poincare-type asymptotic expansions. 

The expansion in (|2.9[) has several advantages compared with (|2.ip. 



In the recursion relation (2.5), both u s and u' s occur in the left-hand side. By using 
computer algebra it is not difficult to compute any number of coefficients u s , but the 
relation for the polynomials 4> s (r) is simpler than (|2.5| ), with this respect. 



The quantity h(n) in (|2.9j ) is defined as an exact relation, and not, as g{n) in (|24|), by 
an asymptotic expansion (cf. (|2.6[)). 
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• Most important, the expansion in (2.9) has a double asymptotic property: it holds if 
one or both parameters t and [i are large, and not only if [i is large. 

For the function V(a, z) we have 

v(-\f,ytj2)= f* P„(t), P M (t) ~ E(-l)^, (2.14) 



where the (j> s {j) are the same as in (|2.9| ). This expansion is a modification of 0-(11.19) (see 
also 0-(2.12)). 

For the derivatives we can use the identities 



(t — 1) 4 ^0 ^ 

Idr^W = V(t a -l)W*Q /l (t), Q„(t) ~ 

ar (< - 1) 4 5=0 ^ 

The coefficients ^ s can be obtained from the relation 



(2.15) 



= Mr) + 2r(r + l)(2r + l)0 s _i(r) + 8r 2 (r + 1) J ^^, (2.16) 



s = 0, 1, 2, . . . . The first few are 



^i(t) = ^(28r 2 + 42r + 15),i; 
r 2 

^(t) =-2gg(7280r 4 + 21840r 3 + 23028r 2 +9684r + 1215), (2.17) 
r 3 

^g(t) = -__(30430400t 6 + 136936800r 5 + 246708000r 4 + 
51o4U 

224494200r 3 + 106122312r 2 + 23489190r + 1658475). 
This gives the modifications (see 0-(4.13)) 

U'{-\^,^ =-JLh(»)(t 2 -l)k-^G,(t), G ,(t)~E^ (2.18) 

and 

Remark 2.1. The functions P«(i), (?«(£), Pu(f) and Q^{t) introduced in the asymptotic 
representations satisfy the following exact relation: 

F^t) Quit) + G^t) P^t) = 2. (2.20) 

This follows from the Wronskian relation (1.6). The relation in ( |2.20| ) provides a conve- 
nient possibility for checking the accuracy in numerical algorithms that use the asymptotic 
expansions of F^it) , G ^(i) , P^t) and Q^(t). 
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2.2 The case a < 0, z < -2y/^a, -a - z > 0. 



For this case we mention the modification of 0-(11.16). That is, for t > 1 + e we have the 
representations 



U (-^ 2 ,-fitV2) 



r(| + ^fj 2 ) cos ^|vr^ 2 
\i v / 7r/i 2 ( / u) 



(2.21) 



where P^(i) and P^(t) have the expansions given in (|2.9|) and (|2.14| ), respectively. An ex- 
pansion for V{—\n 2 ,—^ty/2) follows from the second line in (1.7), ( p.9| ) and (2.21). A few 
manipulations give 



v(-^ 2 ,-fxtV2 



(t 2 - 1)^(1 + 1^) 



r(i + i^)sin Utt/x' 
cos (*„„») ^ F„<t) ^ P„( t ) 



(2.22) 



Expansions for the derivatives follow from the identities in (2.15). If a = — ^fj, 2 = —n — \-,n 



0, 1, 2, . . ., the cosine in ( 2.21 ) vanishes, and, hence, the dominant part vanishes. This is the 
Hermite case, cf. (1.5). 



2.3 The case a < 0, -2\f ZI a < z < 2y/ II a. 

For negative a and — 1 < t < 1 the expansions are essentially different, because now oscillations 
with respect to t occur. We have (0-(5.11) and 0-(5.23)) 



u(-^ 2 ,fitV2 



2<?(a0 



. 2 1 \ (-l) S U2s(t) . ( 2 1 \ST^ (-l) S U 2s+ l(t) 



=0 



s=0 



1-i 



2\3s+§ ,,4s+2 



(2.23) 



with u s (t) defined in (|2.5|) and #(//) in ( |2.6|) , and 

U'(-y 2 ,^tV2) ^ ^9(^(1 -t 2 f x 

(-1)^(0 , / , , \^ (-l) s <; 2s+1 (t) 



sm 



s=0 lJ- ; m s=ol i 



t 2 )3s +§^4 S +2 



(2.24) 



as ^ — > oo, uniformly with respect to \t\ < 1 — e, where the coefficients v s are given by (see 
0-(4.15)) 

v 8 (t) = u s (t) + ±tu a _i(i) - r a _ 2 (t), (2.25) 



and 



?7 = | arccos i — |iV 1 — t 2 . 



(2.26) 
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For the function V(a,z) we have (O-(11.20) and 0(2.12)) 



r(J + J/i 2 )(i-^ 



cos (^r? + Jvrj - sm ^ + g 

l s 

2\i 



s=tr 



t/-/ / i„2 „ + n?\ »V2g{ii)(l-t 



(2.27) 



sin 



r(i + ^ : 

N1 (-l) S V2s(t) 



X 



s=0 ^ ' " s=OV^ 



(-l)^ 2B +l(t) 
2\3s+| ,,4s+2 



By using the Wronskian relation (1.6) it follows that we have the following asymptotic 
identity 

y 1 !.)>.(/) y (-!)>(*) 1 
^(i-^) 3 V 4s s(i-^) a V 4 



s=0 



^ (-l) 8 U 2a +l(t) ~ (-l)^ 2s+1 (t) 



0(1 -t Z ) 



2",,4s+2 



•=0(1 - t 2 )" 



3H 



2 ,.4s+2 



(2.28) 



^2 1 2' 

rvj ■— — 75 

2/X V / 7T5( (//) 



1 



+ 



2021 



576ii 4 2488320// 



2.3.1 Modified expansions 

We can give modified versions based on our earlier modifications, with g(p) replaced with 
h(fi), and so on. Because in the present case t belongs to a finite domain, the modified 
expansions don't have the double asymptotic property. We prefer Olver's versions for this 
case. 

This completes the description of U(a, z), U'(a, z), V(a, z), V'(a, z) in terms of elementary 
functions for negative values of a. 

2.4 The case a > 0, z > 0, a + z > 0. 

For positive values of a the asymptotic behavior is rather simple because no oscillations occur 
now. The results follow from Olver's expansions O-(ll.lO) and 0-(11.12). The modified forms 
are _ 



^-rF,(t), i^)~£(-i) 
{t 2 + 1)4 



s=0 



2s 



where 



£= \ [Wl + t 2 + In (t + Vl + t 2 

1,,2 1..2 1 1 ,,2 1 

The coefficients (j) s in ( |2.29| ) are the same as in (|2.9| ), with r replaced by 



(2.29) 

(2.30) 
(2.31) 



.VTTT 2 



(2.32) 
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For the derivative we have 



tf'(i M 2 ,W2) =~tih(jj l )(l+t i )K-^G a (z), G a (z)~p(-iytjj?l, (2.33) 



where V's( T ) is given in ( p. 16 ), with f defined in ( 2.32| ). 



2.5 The case a > 0, z < 0, a - z > 0. 

Olver's expansion O-(ll.lO) and 0(11.12) cover both cases z > and z < 0. We have the 
modified expansions 



[/' (^ 2 , - = -^^fr^MM)^! + t 2 ) 3 Qfj,(t), 



(2.34) 



where 



In Section 4.1.2 we give details on the derivation of these expansions. 

Remark 2.2. By using the second relation in (1.7), the representations for V(a,z) and 
V'(a, z) for positive a can be obtained from the results for U(a,z) and U'(a,z) in Q2.29Q , 
(gljD and (Ijl). 



Remark 2.3. The functions F^{t) , G M (i) , P M (t) and Q„(i) introduced in ( ggf ), (|3|) and 
( 2.34| ) satisfy the following exact relation: 



F^(t)Q^t) + Gp(t)Pp{t) = 2. (2.35) 
This follows from the Wronskian relation 

U(a, z)U'(a, -z) + U'(a, z)U(a, -z) = - \ 

See also Remark 2.1. 

Remark 2.4. The expansions of Sections 2.4 and 2.5 have the double asymptotic property: 
they are valid if the a + \z\ — > oo. In Sections 2.4 and 2.5 we consider the cases z > and 
z < 0, respectively, as two separate cases. Olver's corresponding expansions O-(ll.lO) and 
0-(11.12) cover both cases and are valid for — oo < t < oo. As always, in Olver's expansions 
large values of \i are needed, whatever the size of t. 

In Figure 1 we show the domains in the t, a— plane where the various expansions of U(a, z) 
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of this section are valid. 




t 



(3.16) (3.1) 
0-(9.7) 0-(8.11) 

Figure 1. Regions for the modified asymptotic expansions of U(a,z) given in 
Section 2 and the Airy-type expansions of Section 3 (which are valid in much 
larger domains than those indicated by the arrows). 



3 Expansions in terms of Airy functions 

The Airy-type expansions are needed if z runs through an interval containing one of the 
turning points ±2a/^o, that is, t = ±1. 



3.1 The case a < 0, z > 0. 

We summarize the basic results 0-(8.11), 0-(8.15) and 0-(11.22) (see also 0-(2.12)): 



Ai(^C) ^(0 + ^^(0 



Hi 



v(-U 2 ,^tV2) = 



<KC) 



Ai(/iaC) 



4 

M 3 



c M (C) + AiVC)^(C) 



2^3ff(/i)0(C) 

r(i + V) 

1 2 

(2tt) /i3g(/i) 

U 1„2^ 



Bi(M*0 ^(0 + 5^ B„(C) 



J 0(C)r(i + i/i 



Bi^K) 



4 

Ai3 



(3.1) 

(3.2) 
(3.3) 
(3.4) 



The coefficient functions A^Q, B^Q, C M (£) and -D M (C) have the following asymptotic ex- 
pansions 

s=0 ^ s=0 ^ 
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cwo~E#, ^(O-E^ (3-6) 

as fj, — > oo, uniformly with respect to t > — 1 + 6, where 5 is a small fixed positive number. 
The quantity ( is defined by 

I(-C) § =v(t), -Kt<l, (C<o) 

§c^ =m, i<t, cc>o) 

where 77, £ follow from (2.26), (2.2), respectively, and 

^(0=(^) i - (3.8) 

The function £(t) is real for t > —1 and analytic at t = 1. We can invert into t(C) ; an d 
obtain 

t = 1 + 2~V3 C - ^2" 2 / 3 C 2 + ^C 3 + • • • • 

The function g(fj,) has the expansion given in (2.6) and the coefficients a s ((),b s (() are given 
by 

2s _3 2s+l a 

a s (Q = J2 PmC ^ m A 2 s- m {t) VCbs(C) = ~J2 a mC i * m A 2s -m+l(t), (3.9) 
m=0 m=0 

where A s (t) are used in (2.1), cxq = 1 and 

_ (2m + l)(2m + 3)---(6m-l) 6m + 1 

a -" m!(144)- ' 6m-l° m - ^ Uj 

A recursion for a m reads: 

(6m + 5)(6m + 3)(6m+ 1) 
am+1= ° m 144(m + l)(2m + l) ' - = ' 1 ' 2 '-- 

The numbers a m ,f3 m occur in the asymptotic expansions of the Airy functions, and the 
relations in ( |3.9j ) follow from solving ( |3.1[) and (|3.3| ) for Au(C) an d £?„(£), expanding the Airy 
functions (assuming that £ is bounded away from 0) and by using (2.1) and a similar result 
for V(a,z) (0-(11.16) and 0-(2.12)). 

For negative values of C (i.e., — 1 < t < 1) we can use (0-(13.4)) 



O»(0 =(-l) S E™=0/3 m (-C) 2 ^2s- m (t), 

C& S (C) = (-l)*" 1 E^o 1 « m (-C)" §m ^2s- m+ i(t), 



(3.11) 



where 



(i-t 2 )t s 

The functions C M (C) and D^(() of ( |3^ ) and ( |3.4| ) are given by 



c„(C) = x(C)4(0 + ^(C) + (B^O, D^O = 4(0 + 4 fxCCWO + s;(c)l . (3.12) 
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The coefficients c s (C) and d s (() in ( |3.6| ) are given by 

c s (0 = x(CMC) + o'.(C) + CMC). <*«(C) = o.(C) + x(C)&-i(C) + d(C), (3-13) 

where 

0'(C) i - 2#(Q] 6 

* (c) = 0(0 = ic ' (3 - 14) 

with c/)(C) given in (|3.8|). Explicitly, 



, 2s+l q 2s 

1 x — ^ — Km ^ — ^ — xm 



/3 m C #2 S -m+l(r) 4(C) = - 2^ «mC ^2s-w(t), (3-15) 
v ^ m=0 m=0 

3 S 

where S s (r) = v s (t)/(t 2 — 1)^ , with v s (t) defined in (2.25). Other versions of (|3.15| ) are 
needed for negative values of C4-6-, if — 1 < £ < 1; see (|3.11| ). 

3.2 The case a < 0, z < 0. 

Near the other turning point t = —1 we can use the representations (0-(9.7)) 



U (-^ 2 ,-fitV2) =2^^g{ii)<t>{Q sin(j7r//) Ud0QA^Q+ 

^J^^(C)} + cos (iV) + W 



(3.16) 



as \i — ► oo, uniformly with respect to t > —1 + 5, where <5 is a small fixed positive number. 
Expansions for V(a,z) follow from ( |3.1[ ) and ( 3.16[ ) and the second relation in (1.7). Results 



for the derivatives of U(a, z) and V(a, z) follow easily from the earlier results. 
3.3 Modified forms of Olver's Airy-type expansions 



Modified versions of the Airy-type expansions fl3.1| ) - Q3.4| ) can also be given. In the case of 
the expansions in terms of elementary functions our main motivation for introducing modi- 
fied expansions was the double asymptotic property of these expansions. In the case of the 
Airy-type expansions the interesting domains for the parameter t, from a numerical point 
of view, are finite domains that contain the turning points ±1. So, considering the expan- 
sions given so far, there is no need to have Airy- type expansions with the double asymptotic 
property; if \x remains finite and \t\ S> 1 we can use the expansions in terms of elementary 
functions. However, we have another interest in modified expansions in the case of Airy-type 
expansions. We explain this by first discussing a few properties of the coefficient functions 
4,(0,^(0,(^(0 and £>„(C). 

By using the Wronskian relation (1.6) we can verify the relation 

i m + iv 2 ) 

.4,(0^(0 - ^«)C„(C) = ^j^-y (3.17) 

where g{n) is defined by means of an asymptotic expansion given in (2.6). By using the 
differential equation (0-(7.2)) 

^=[ M 4C + *(C)]W, (3.18) 
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where 

*(C) 



5_ _ (3t 2 + 2)( = 2l/3 



16C 2 4(t 2 - l) 3 



_9_ , _J_ 2 -l/3f 
280 ~ l-W ^ 



ir»o 



1359 q-2/3a2 i 196 a3 
26950 ^ ^ 8125^ 



we can derive the following system of equations for the functions A^((), 

A" + 2C-B' + B - V(()A = 0, 
B" + 2^ A' -*(C)B =0, 



(3.19) 



(3.20) 



where primes denote differentiation with respect to £. A Wronskian for this system follows 
by eliminating the terms with ^(C)- This gives 



2/x 4 A 1 'A + AB" - A"B - 2(B'B - B 2 



0, 



which can be integrated: 



1 + V) 



M 4 4(C) + ^(CK(C) - A'JQMO - mo = ^ 



(3.21) 



2jj,^/TTg 2 {^y 

where the quantity on the right-hand side follows from ( |3.17| ) and ( p.!2| ) . It has the expansion 



1 



+ 



2021 



576^ 4 2488320/z 8 



(3.22) 



as follows from 0-(2.22) and 0-(5.21). 

As mentioned before, the interesting domain of the Airy-type expansions given in this 
section is the domain that contains the turning point t = 1, or £ = 0. The representations of 
the coefficients of the expansions given in ( [3.9D cannot be used in numerical algorithms when 
|£| is small, unless we expand all relevant coefficients in powers of £. This is one way how to 
handle this problem numerically; see fl8|| . In that paper we have discussed another method 
that is based on a system like ( 3,2C| ), with applications to Bessel functions. In that method the 
functions A^(() and B^Q are expanded in powers of (, for sufficiently small values of |£|, say 
|C| < 1, and the Maclaurin coefficients are computed from ( 3.20| ) by recursion. A normalizing 
relation (the analogue of ( 3.21 )) plays a crucial role in that algorithm. The method works 
quite well for relatively small values of a parameter (the order of the Bessel functions) that 
is the analogue of 

When we want to use this algorithm for the present case only large values of ji are allowed 



because the function g(n) that is used in (|3.l|) - (|3.4|) and (|3.2l| ) is only defined for large values 
of [i. For this reason we give the modified versions of Olver's Airy-type expansions. The 
modified versions are more complicated than the Olver's expansions, because the analogues of 
the series in ( |3.5[) and ( |3.6[) are in powers of ^~ 2 , and not in powers of n~ A . Hence, when we use 
these series for numerical computations we need more coefficients in the modified expansions, 
which is certainly not desirable from a numerical point of view, given the complexity of 
the coefficients in Airy-type expansions. However, in the algorithm based on Maclaurin 
expansions of the analogues of the coefficient functions A^((), B^((), C^(() and D^(() this 
point is of minor concern. 

The modified expansions are the following: 



U 



-i/i 2 ,/xt\/2 



rg + V)<MC) 

li*h(n) 



AfVK) 



(3.23) 
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V (-^^tV^j 



0(0 



/X3 



The functions -^(C) an d G>(C) have the following asymptotic expansions 



(3.24) 



^(0~E#' G m(0~E#- (3-25) 

s=0 V s=0 V 

The quantity £ and the functions 0(C) an d 6,(a*) are as in Section 3.1. Comparing ( 3.23j ), 
(|3~2l ) with Q, (H) we conclude that 



F M (C)=F(^(C), G M (C) = H^B^C), %) = ^W " ^ 

1 ^2 ' 2^ 



The function #"(//) can be expanded (see 0-(2.22), 0-(2.27), 0-(6.2) and (2.6)) 

oo 

^)~l + IE(-l) S 7Tl)I 5 (3-27) 
where 7 S are the coefficients in the gamma function expansions 



The first few coefficients are 

( 1 1 1003 

7o - , 7i - -^4, 72 - Y?^' ^ 3 ~ 414720 ' 

The second expansion in ( |3.28| ) can be used in ( 0.26 ) to find relations between the coeffi- 
cients a s (C) and 6 S (C) °f (|3.5|) and of / S (C) and </ s (C) of ( 3.25[ ). That is, 



/o(C) = i, /i(C) = 4 /2(C) = 01(0 + ^, M0 = i-MO 



1003 



24 1 J^vs/ J- vs / 5T6 ' 24 103680' 

90(C) = 60(C), 91(C) = o%(0, 92(C) = 61(C) + ^feo(C), 93(C) = ^6i(C) - tStMO- 



The coefficients / s (C)>ffs(C) ca n also be expressed in terms of the coefficients s (r) that are 
introduced in (2.9) by deriving the analogues of (|3.9|). 
The system of equations ( 3.20| ) remains the same: 



F" + 2CG' + G - *(C)F =0, 
G // + 2/F / -^(C)G =0, 



(3.29) 



and the Wronskian relation becomes 



A?(C) + ^(0^(0 - ^(C)^(C) - CG?(C) = / ^ff^j - (3.30) 

The right-hand side has the expansion (see ( 3.28Q and (2.7)) /x 4 J2^=o(~^Yls/ {\^ 2 ) s ■ Observe 
that ((IJTJ) is an exact relation, whereas ( |3.2l| ) contains the function g(fi), of which only an 
asymptotic expansion is available. 
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3.4 Numerical aspects of the Airy-type expansions 

In |18| , (Section 4), we solved the system (|3.29|) (for the case of Bessel functions) by substi- 
tuting Maclaurin series of F(Q, G(C) and ^ (C)- That is, we wrote 

oo oo oo 

n=0 n=0 n=0 



where the coefficients i[) n can be considered as known (see (|3.19| )), and we substituted the 
expansions in ( |3.29| ), This gives for n = 0, 1, 2, . . . the recursion relations 



(n + 2)(n + l)c n+2 + (2n + l)d n = p n , p n = ^V'fcCn-fc, 



k=0 



(3.31) 



(n + 2)(n + l)d n+2 + 2/j, 4 (n + l)c„+i = a n , a n = ^ip k d n _ k . 

k=0 



If /x is large, the recursion relations cannot be solved in forward direction, because of numerical 
instabilities. For the Bessel function case we have shown that we can solve the system by 



iteration and backward recursion. The relation in ( |3.30| ) can be used for normalization of the 
coefficients in the backward recursion scheme. 

For details we refer to [jig]. The present case is identical to the case of the Bessel functions; 
only the function ^(C) is different, and instead of p? in (3.31) we had the order v of the Bessel 
functions. 

4 Expansions from integral representations 

The expansions developed by Olver, of which some are given in the previous sections, are 
all valid if \a\ is large. For several cases we gave modified expansions that hold if at least 
one of the two parameters a, z is large and we have indicated the relations between Olver's 
expansions and the new expansions. The modified expansions have in fact a double asymp- 
totic property. Initially, we derived these expansions by using integral representations of the 
parabolic cylinder functions, and later we found the relations with Olver's expansions. In this 
section we explain how some of the modified expansions can be obtained from the integrals 
that define U(a, z) and V(a, z). Again we only consider real values of the parameters. 

4.1 Expansions in terms of elementary functions by using integrals 

4.1.1 The case a>0,z>0; a + z > 0. 

We start with the well-known integral representation 

e f°° a— A — Au> 2 — zw , 

U(a,z) = — / w 2 e dw, a > -± (4.1) 

T(a + Jo z 

which we write in the form 

/+i "i* 3 ,oo _i , 

U(a, z) = - / w 2 e- z *W dw, (4.2) 

T(a + ^) Jo 
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where 



(w) = w + \w 2 - Ainu;, A=-^. (4.3) 



The positive saddle point wq of the integrand in (4.3) is computed from 

d(j)(w) w 2 + w — A 



dw w 

giving 



0, (4.4) 



w = | 



Vl + 4A - 1 . (4.5) 



We consider z as the large parameter. When A is bounded away from we can use Laplace's 
method (see [|ll]] or |^2| ). When a and z are such that A —* Laplace's method cannot be 
applied. However, we can use a method given in [15] that allows small values of A. 

To obtain a standard form for this Laplace-type integral, we transform w — > t (see fl6f ) 
by writing 

=t- Alnt + A (4.6) 

where A does not depend on t or u;, and we prescribe that w = should correspond with 
i = and w = wq with i = A, the saddle point in the t— plane. 
This gives 

a+l -Xz 2 -Az 2 

U(a,z) = — — / t ^e~ zt f{t)dt, (4.7) 



(l + 4A)4r(a+^)^o 
where 

/(*) = (1 + 4A) 4 #^ = (1+4A)*# 2 t " A r - (4.8) 
' /wv ' \l w dt y J \l t w 2 + w- X v/ 

l 

By normalizing with the quantity (1 + 4A) ? we obtain /(A) = 1, as can be verified from ( [4 . 8[) 
and a limiting process (using l'Hopital's rule). The quantity A is given by 

A = \wl + w - Alnw - A + A In A. (4.9) 
A first uniform expansion can be obtained by writing 

CO 

/(t) = 5> n (A)(i-Ay\ (4.10) 

n=0 

Details on the computation of a n (A) will be given in the appendix. 
By substituting (|4.10| ) into (0) we obtain 

U(a, z) ~ — ^ V a n (A)P„(a)z- 2 ", (4.11) 



*"" 2 (l + 4A)i n=0 
where 

_2a+2n+l /-co „_ 1 

Z I 4. 2 „-2 2 t/ 



P n (a) = —— T - t 2 e- zt (t-X) n dt, n = 0,1,2,... . (4.12) 
r(a + Jo 

The P n (a) are polynomials in a. They follow the recursion relation 

P n +i(a) = (n + |)P n (a) + anP n _i(a), n = 0, 1,2, . . . , 
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with initial values 

P (a) = 1, Pi(a) = \. 
We can obtain a second expansion 



U(a,z) 



a+4 



(4.13) 



;i + 4A)3 fc=o 



with the property that in the series the parameters A and z are separated, by introducing a 
sequence of functions {/&} with fo(t) = f(t) and by defining 



Vt 



fk( t)-fk(X) 
t-X 



, fe = 0,1,2, 



(4.14) 



The coefficients fk{X) can be expressed in terms of the coefficients a n (X) defined in ( 4.10| ). 
To verify this, we write 



A(t) = £ag' ) (A)(*-A) B . 

n=0 

and by substituting this in ( 4.14| ) it follows that 

ai k+1 \X) = X(n + 1)4% (A) + (n + |) 4% (A), fc > 0, 

Hence, the coefficients /fc(A) of ( 4.13j ) are given by 

/ fc (A) = 4 fc) (A), k>0. 



n > 0. 



We have 



/o(A) =1, 

/i(A) =i[a 1 (A) + 2Ao 2 (A)], 

/ 2 (A) = | [l2A 2 a 4 (A) + 14Aa 3 (A) + 3a 2 (A)] , 

/ 3 (A) = | [l20A 3 a 6 (A) + 220A 2 a 5 (A) + 116Aa 4 (A) + 15o 3 (A)] . 



Explicitly, 



/o(A) 
/i(A) 

/ 2 (A) 

/ 3 (A) 



24 

2 



(20a 2 - 10a - 1), 



1152 



(6160a 4 - 6160a 3 + 924a 2 + 20a + 1), 



P 



-(27227200a 6 - 40840800a 5 + 16336320a 4 - 



where 



a 



1 + 



414720 

1315160a 3 - 8112a 2 + 2874a + 1003), 



(2a - l) 2 2z 2 



\/4a + z 1 



a V4a + z 2 (z + V4a + z 2 ) ' 



We observe that /fc(A) is a polynomial of degree 2k in a multiplied with p k . 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



(4.19) 



(4.20) 
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If a and z are positive then a E [0, 1]. Furthermore, the sequence {p k / z 2k } is an asymptotic 
scale when one or both parameters a and z are large. The expansion in ( |4.13] ) is valid for 
z — > oo and holds uniformly for a > 0. It has a double asymptotic property in the sense that 
it is also valid as a — > oo, uniformly with respect to z > 0. As follows from the coefficients 
given in (4.19) and relations to be given later, we can indeed let z — ► in the expansion. 

The expansion in fl4.13| ) can be obtained by using an integration by parts procedure. We 
give a few steps in this method. Consider the integral 



r(a + ^)7o 



-zH 



f(t) dt, 



(4.21) 



We have (with A = a/z 2 
/(A) 



FJz) 



Via 



ttJo f 



1 ™ - a ~h„-z 2 t 



1 i> 

-2a-l 



r(a + i) 



[f(t)-f(X)}dt 



/(A) 



-2a-l 



1 fCXD , 2" 

32~F7^ i 1 WO r 



i [/(*) - /(A)] 



z 2 r(a + 
/(A) + 



t- A 



de -z 2 (t-Xlnt) 



1 



^r(a + 



-zH 



fi(t)dt, 



where f\ is given in ( [4.14 ) with /o = /. Repeating this procedure we obtain ( 4.1 3| ) . More 
details on this method and proofs of the asymptotic nature of the expansions (4.11) and Q4.13 ) 
can be found in our earlier papers. We concentrate on expansion ( 4.13 ) because ( 4.11 ) cannot 
be compared with Olver's expansions. 

To compare ([OJD with Olver's expansion (2.16), we write 



\i\fit. 



(4.22) 



Then the parameters a and p defined in fl4.20| ) become 

t 



a 



1 + 



T+l, P 



2t 2 



vT+t^t + VT+t 2 ') ' 



(4.23) 



where r is given in (2.32). After a few manipulations we write ([OJD in the form (cf.(2.29)) 



Uiy 2 ,iity/2 



(t 2 + 1) 3 



k=0 



2k 



where 



and 



i = \ [Wl + t 2 + In (t + Vl + t 2 ) 



1,,2 1 ,,2 1 1„2 1 

h(p) = e 4 n 2 2 2 4 4 



(-l) fc 



(4.24) 

(4.25) 
(4.26) 

(4.27) 
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Explicitly, 



01 (a) = ^^(20a 2 -10(7-1), 

^((t) = (6160a 4 - 6160a 3 + 924a 2 + 20a + 1), (4.28) 

zoo 

03(°") = ^!-~^ (27227200a 6 - 40840800a 5 + 16336320a 4 - 
51o4U 

1315160a 3 - 8112a 2 + 2874a + 1003), 

where a is given in Q4.23D . Comparing ( |4.24| ) with (2.29) we obtain 0fc(a) = <pk(T),k > 0, 
because a = 1 + f. 

4.1.2 The case a > 0, z < 0; a - z > 0. 

To derive the first expansion in (2.34) we use the contour integral 

A/27TP 4 " r(a+A) /■ zs +ls 2 -a-A 

r(a + £) 27T? Jc 

where C is a vertical line in the half plane 5Rs > 0. This integral can be transformed into a 
standard form that involves the same mapping as in the previous subsection. We first write 
(by transforming via s = zw) 

. . z a r(a + i) r z 2 { w +W) - a - A , « a r(a + ^) f , dw 

H a (z) = — / e 1 w 1 dw = — / e z W°> (4.30) 

2m Jr 2m Jr \/W 



where <f>(w) is defined in ( [4.3[ ). By using the transformation given in ( fig ) it follows that 

z^ a T{a + \)e Az2 , 2 _ Q _l 

H a (z) = — / e z \ 1 f(t) dt. (4.31) 

2m Jc 



The integration by parts method used for ( [4.21 ) gives the expansion (see 



a Az 2 oo f,(\\ 

^)~7ff^rE(-l) fc ^' (4-32) 
(4a + z 2 )4 fc=0 z 



where the /fc(A) are the same as in fl4.13|) ; see also ( 4.18 ). This gives the first expansion of 
(2.34). 



Remark 4.1. The first result in (2.34) can also be obtained by using (|4.l| ) with z < 0. The 
integral for U[a, —z) can be written as in ( |4.2| ), now with cj)(w) = \w 2 — w — In A, A = a/z 2 . 

In this case the relevant saddle point at Wq = (1 + y/1 + 4A)/2 is always inside the interval 
[1, oo) and the standard method of Laplace can be used. The same expansion will be obtained 
with the same structure and coefficients as in (2.34), because of the unicity of Poincare-type 
asymptotic expansions. See also Section 4.1.4 where Laplace's method will be used for an 
integral that defines V(a,z). 
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4.1.3 The case a<0,z>2y/^a, -a + z>0. 



Olver's starting point (2.1) can also be obtained from an integral. Observe that ( |4.1[) is not 
valid for a < —\. We take as integral (see [0], p. 687, 19.5.1) 



U(-a,z) 



r(i + a) l : 



2-ni 



.1.2 



(4.33) 



where a is a contour that encircles the negative s— axis in positive direction. Using a trans- 
formation we can write this in the form (cf. (|4.2|) ) 



U(-a,z) 



T(i + a) i a -W 

z l e 4 

2m 



w dw, 



where 



(w) = W — t;W 2 — Ainu;, A = 



The relevant saddle point is now given by 

w = | (l - Vl - 4A) , < A < \. 



(4.34) 
(4.35) 

(4.36) 



When A — > the standard saddle point method is not applicable, and we can again use 
the methods of our earlier papers (]l5[| and |l6[) and transform 

(4.37) 



(w) = t- \lnt + A, 



where the points at — oo in the w— and t— plane should correspond, and w = w$ with t = A. 
We obtain 



U(-a,z) 



r(i + o) 1 a -l z z+z 2 A 

-z z e 4 



1 



f(t)dt, 



(l-4A)4 2vri 

where a is a contour that encircles the negative t— axis in positive direction and 

l 



/(i) = (l-4A) 



t- A 



1 T dw 

w dt \ t w — w L — X 



(4.38) 



(4.39) 



Expanding f(t) as in ( [4.10| ), and computing /fc(A) as in the procedure that yields the relations 
in fl4.18p , we find that the same values /fc(A) as in ( |4.19j ), up to a factor(— l) k and a different 
value of r and p. By using the integration by parts method for contour integrals [15|, that 
earlier produced ( 4.32j ), we obtain the result 



U(-a,z) 



a Az2 ~b 2 oo 

z a e 



r EC" 1 )' 



(z 2 -4a)4 fc=0 



where the first /^(A) are given in ( 4.19 ) with 

- (2"-l) a 



(7 



1 + 



aZ-z 2 - 4a 



2z 2 



(4.40) 



V^ 2 - 4:a+ (z + Vz 2 - 4a) ' 



(4.41) 



This expansion can be written in the form (2.9). 
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4.1.4 The case a < 0, z < -2> 

We use the relation (see (1.7)) 



-a - z > 0. 



U (—a, —z) = sin7ra U (—a, z) + 



7T 



r(l 



V(-a,z), 



(4.42) 



and use the result of U(—a,z) given in (|4.40| ) or the form (2.9). An expansion for V(—a,z) 
in Q4.42p can be obtained from the integral (see Q) 



1,2 

3 



V(a,z) 



2tt 



-i- 9 



s +2s a- 
s 



(is, 



(4.43) 



71U72 



where 71 and 72 are two horizontal lines, 71 in the upper half plane 9s > and 72 in the 
lower half plane 9s < 0; the integration is from 5fts = — 00 to !Rs = +00. (Observe that when 
we integrate on 71 in the other direction (from $ts = +00 to $ts = 00) the contour 71 U 72 
can be deformed into a of ( [4,33j ), and the integral defines U(a,z), up to a factor.) We can 
apply Laplace's method to obtain the expansion given in (2.14) (see Remark 4.1). 



4.2 The singular points of the mapping ( f4.6|) 

The mapping defined in (4.6) is singular at the saddle point 



W- 



1 + 4A + 1 



(4.44) 



If A = then w. 
the estimate: 



-1 and the corresponding t— value is —\. For large values of A we have 

0.4356" 



t(w-) ~ A 



-0.2785 



(4.45) 



This estimate is obtained as follows. The value t— = t(v)—) is implicitly defined by equation 
flQjl) with w = W-. This gives 



t- - Alni_ - A + A In A 



Wl + 4A± Avri + A In 



4A 



±\iri - 2y/\ 



1 



(1 + TTTlA) 5 



1+ 24A + °^ 



(4.46) 



as A — > 00. The numerical solution of the equation s — Ins — 1 = ±ni is given by s± = 
0.2785- • • e^ m . This gives the leading term in (4.46). The other term follows by a further 
simple step. 



4.3 Expansions in terms of Airy functions 

All results for the modified Airy-type expansions given in Section 3.3 can be obtained by using 
certain loop integrals. The integrals in (4.33) and (4.43) can be used for obtaining (3.23) and 
(3.24), respectively. The method is based on replacing 4>(w) in (4.34) by a cubic polynomial, 
in order to take into account the influence of both saddle points of <f>(w). This method is first 
described in ||; see also [0] and 
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5 Numerical verifications 



We verify several asymptotic expansions by computing the error in the Wronskian relation 
for the series in the asymptotic expansions. Consider Olver's expansions of Section 2.3 for 
the oscillatory region — 1 < t < 1 with negative a. We verify the relation in (2.28). Denote 
the left-hand side of the first line in (2.28) by W(fJ,, t). Then we define as the error in the 
expansions 



1 



+ 



2021 



576^ 4 2488320/x fc 



(5.1) 



Taking 3 terms in the series of (2.23), (2.24) and (2.27), we obtain for several values of \x 
and t the results given in Table 5.1. We clearly see the loss of accuracy when t is close to 1. 
Exactly the same results are obtained for negative values of t in this interval. 





5 


10 


25 


50 


100 


t 

.00 


.32e-09 


.78e-13 


.13e- 


17 


.32e- 


21 


.78e-25 


.10 


.26e-09 


.63e-13 


.lle- 


17 


.26e- 


21 


.63e-25 


.20 


.81e-10 


.20e-13 


.33e- 


18 


.82e- 


22 


.20e-25 


.30 


.16e-08 


.39e-12 


.65e- 


17 


.16e- 


20 


.39e-24 


.40 


.88e-08 


.22e-ll 


.36e- 


16 


.89e- 


20 


.22e-23 


.50 


.51e-07 


.13e-10 


.21e- 


15 


.52e- 


19 


.13e-22 


.60 


.40e-06 


.99e-10 


.17e- 


14 


.40e- 


18 


.99e-22 


.70 


.53e-05 


.13e-08 


.22e- 


13 


.54e- 


17 


.13e-20 


.80 


.20e-03 


.50e-07 


.84e- 


12 


.20e- 


15 


.50e-19 


.90 


.35e-00 


.24e-04 


.41e- 


09 


.10e- 


12 


.25e-16 



Table 5.1. Relative accuracy A(/i,i) defined in (|5.1| ) for the asymptotic series 
of Section 2.3. 



Next we consider the modified expansions of Section 2.1. Denote the left-hand side of 
(2.20) by W(fi,t). Then we define as the error in the expansions 



A(/i,t) :-- 



(5.2) 



When we use the series in (2.9), (2.14), (2.18) and (2.19) with 5 terms, we obtain the results 
given in Table 5.2. We observe that the accuracy improves as /j, or t increase. This shows the 
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double asymptotic poperty of the modified expansions of Section 2.1. 



10 



25 



50 



100 



t 

1.1 


.51e-01 


.48e-05 


.72e- 


10 


.18e- 


13 


.43e-17 


1.2 


.39e-04 


.79e-08 


.13e- 


12 


.32e- 


16 


.78e-20 


1.3 


.83e-06 


.19e-09 


.32e- 


14 


.78e- 


18 


.19e-21 


1.4 


.56e-07 


.13e-10 


.23e- 


15 


.55e- 


19 


.13e-22 


1.5 


.71e-08 


.17e-ll 


.29e- 


16 


.70e- 


20 


.17e-23 


2.0 


.10e-10 


.25e-14 


.43e- 


19 


.10e- 


22 


.25e-26 


2.5 


.21e-12 


.52e-16 


.87e- 


21 


.21e- 


24 


.52e-28 


5.0 


.12e-16 


.28e-20 


.48e- 


25 


.12e- 


28 


.28e-32 


10.0 


.20e-20 


.48e-24 


.81e- 


29 


.20e- 


32 


.48e-36 


25.0 


.30e-25 


.73e-29 


.12e- 


33 


.30e- 


37 


.73e-41 



Table 5.2. Relative accuracy A([i,t) defined in 
of Section 2.1. 



for the asymptotic series 



Finally we consider the expansions of Sections 2.4 and 2.5. Let the left-hand side of (2.35) 
be denoted by W(fi,t). Then we define as the error in the expansions 



AGM) ■-- 



hw(ji,t) 



i 



(5.3) 



When we use the series in (2.29), (2.33) and (2.34) with 5 terms, we obtain the results of 
Table 5.3. We again observe that the accuracy improves as /U or t increase. This shows the 
double asymptotic poperty of the modified expansions of Sections 2.4 and 2.5. 





5 


10 


t 








.00 


.32e-09 


.78e- 


13 


.25 


.12e-09 


.28e- 


13 


.50 


.45e-ll 


.lle- 


14 


.75 


.57e-ll 


.14e- 


14 


1.0 


.27e-ll 


.65e- 


15 


1.5 


.29e-13 


.70e- 


17 


2.0 


.20e-13 


.48e- 


17 


2.5 


.43e-14 


.lle- 


17 


5.0 


.45e-17 


.lle- 


20 


10.0 


.16e-20 


.38e- 
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Table 5.3. Relative accuracy A(//, t) 
of Sections 2.4 and 2.5. 



25 50 100 



13e- 


17 


.32e- 


21 


.78e- 


25 


47e- 


18 


.12e- 


21 


.28e- 


25 


19e- 


19 


.46e- 


23 


.lle- 


26 


24e- 


19 


.58e- 


23 


.14e- 


26 


lle- 


19 


.27e- 


23 


.65e- 


27 


12e- 


21 


.29e- 


25 


.70e- 


29 


81e- 


22 


.20e- 


25 


.48e- 


29 


18e- 


22 


.43e- 


26 


.lle- 


29 


18e- 


25 


.45e- 


29 


.lle- 


32 


64e- 


29 


.16e- 


32 


.38e- 


36 



in ( |5.3| ) for the asymptotic series 



6 Concluding remarks 

As mentioned in Section 1.1 of the Introduction, several sources for numerical algorithms 
for evaluating parabolic cylinder functions are available in the literature, but not so many 
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algorithms make use of asymptotic expansions. The paper |T0[ is a rich source for asymptotic 
expansions, for all combinations of real and complex parameters, where always \a\ has to be 
large. There are no published algorithms that make use of Olver's expansions, although very 
efficient algorithms can be designed by using the variety of these expansions; j| is the only 
reference we found in which Olver's expansions are used for numerical computations. 

We started our efforts in making algorithms for the case of real parameters. We selected 
appropriate expansions from Olver's paper and for some cases we modified Olver's expansions 
in order to get expansions having a double asymptotic property. A serious point is making 
efficient use of the powerful Airy-type expansions that are valid near the turning points of the 
differential equation (and in much larger intervals and domains of the complex plane) . In par- 
ticular, constructing reliable software for all possible combinations of the complex parameters 
a and z is a challengeing problem. 

A point of research interest is also the construction of error bounds for Olver's expansions 
and the modified expansions. Olver's paper is written before he developed the construction of 
bounds for the remainders, which he based on methods for differential equations, and which 
are available now in his book [ p.l|| . 



7 Appendix: Computing the coefficients /fc(A) of (4.13) 

We give the details on the computation of the coefficients fk(X) that are used in (4.13). The 
first step is to obtain coefficients dk in the expansion 

w = d + d 1 (t-X)+d 2 (t-X) 2 + ..., (7.1) 

where do = wq. From (4.6) we obtain 

dw w t — A 



dt t w 2 + w — A 

Substituting fl7.l| ) we obtain 



(7.2) 



^ = wt^v (7 - 3) 

where the saddle point wq is defined in (4.5). From the conditions on the mapping (4.6) 
it follows that d\ > 0. Higher order coefficients dk can be obtained from the first ones by 
recursion. 

When we have determined the coefficients in ( |7.1[ ) we can use (4.8) to obtain the coefficients 
a n (A) of (4.10). 

For computing in this way a set of coefficients fk(X), say /o(A), . . . , /15(A), we need more 
than 35 coefficients d^ in (|7.1| ) . Just taking the square root in (|7.3|) gives for higher coefficients 
dk very complicated expressions, and even by using computer algebra programs, as Maple, 
we need suitable methods in computing the coefficients. 

The computation of the coefficients dk,a n (X) and A (A) is done with a new parameter 
9 £ [0, \tt) which is defined by 

4A = tan 2 9. (7.4) 

We also write 

a = cos 2 ±0, (7.5) 
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which is introduced earlier in (4.20) and (4.23). Then 



l-<r . <t(1-<t) 2a -1 

W0 = 2^rT' A= (2^Tp dl = ^- (7 - 6) 

In particular the expressions for wo and d\ are quite convenient, because we can proceed 
without square roots in the computations. Higher coefficients dk can be obtained by using 

The first relation / (A) = a[, 0) (A) = 1 easily follows from (4.3), (4.8), Q and Q: 



i 
*3 



A 



/o(A) = (1 + 4A) 4 J — di 
V 

Then using (4.8) we obtain 

cos 2 6»(l + 2c) 2 , x . cos 4 6»(20c 4 + 40c 3 + 30c 2 + I2c + 3) 
oo(A) = 1, ox(A) = ~ 6(c + 1)c2 , 02(A) = 24(c + l) 2c 4 • 

where c = \fa = cos \Q. Using the scheme leading to (4.17) one obtains the coefficients /fc(A). 
The first few coefficients are given in (4.19). 

We observe that /fe(A) is a polynomial of degree 2k in a multiplied with p k . If a and z are 
positive then a £ [0, 1]. It follows that the sequence {p k / z 2k } is an asymptotic scale when one 
or both parameters a and z are large, and, hence, that {fk(X)/z 2h } of (4.13) is an asymptotic 
scale when one or both parameters a and z are large. 

Because of the relation in (4.27) and 4>k(c) = 4>k{^)i higher coefficients /fc(A) can also 
be obtained from the recursion relation (2.11), which is obtained by using the differential 
equation of the parabolic cylinder functions. 
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